Strong Three-magnon Scattering in Cuprates by Resonant X-rays 
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We show that Resonant Inelastic X-ray scattering (RIXS) is sensitive to three-magnon excitations 
in cuprates. Even if it requires three electrons to simultaneously flip their spin, the RIXS tri-magnon 
scattering amplitude is not small. At the Cu L-edge its intensity is generally larger than the bi- 
magnon one and at low transferred momentum even larger than the single-magnon intensity. At the 
copper M-edge the situation is yet more extreme: in this case three-magnon scattering is dominating 
over all other magnetic channels. 
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Introduction. Traditionally, inelastic neutron scatter- 
ing is the technique of choice to measure the dispersion of 
elementary magnetic excitations in magnetic materials. 
Only very recently photon scattering became competi- 
tive in this respect, with Resonant Inelastic X-ray Scat- 
tering (RlXSj^ getting energy-sensitive enough to ob- 
serve magnon excitation^^Hll _ experimental endeavor 
underpinned by theorjEHU! RIXS is a photon in - pho- 
ton out technique in which the energy of the incoming 
x-ray photon is tuned to an absorption edge of the ma- 
terial under study. By tuning to the copper L^-edge, for 
instance, the single-magnon dispersion in La2Cu04 was 
recently mapped out and found to be in perfect agree- 
ment with the neutron scattering datcP. At the Cu K- 
edge the RIXS process is indireclfSl and single-magnon 
scattering is forbidden. Instead^ bi-magnon creation is 
allowecplf^ and indeed observed'^. In the latter case mag- 
netic RIXS is thus complementary to optical Raman scat- 
tering, which also measures the bi-magnon, but at zero 
momentum transfer onljEl. 

Here we show that in addition to these rather fa- 
miliar magnetic modes, RIXS also has the unique po- 
tential to examine properties of a virtually unexplored 
magnetic excitation that cannot be probed directly by 
any other technique: the tri-magnon. Tri-magnons are 
three-magnon states that have been identified as im- 
portant and generic features of gapful, low-dimensional 
quantum spin systemJ^. Theoretically it has long been 
known that three reversed spins can form a tri-magnon 
boundstat^i^^Hi^, but the spectroscopic means to directly 
create and probe such an excitation has been lacking so 
far. We find that in RIXS the tri-magnon scattering am- 
plitude is not small, even if it requires three electrons 
to simultaneously flip their spin. At the Cu L-edge the 
three-magnon intensity is generally larger than the bi- 
magnon one and at low transferred momentum q we find 
that it is even larger than single-magnon intensity. At the 
copper M-edge the situation is yet more extreme, as in 
this case the three-magnon scattering channel is the dom- 
inant one. It is thus clear that for a quantitative inter- 
pretation of magnetic L- and Af-edge RIXS data three- 
magnon scattering is indispensable. At the same time the 
presence of a substantial tri-magnon scattering amplitude 



illustrates that RIXS has an unprecedented potential to 
measure properties of interesting multi-particle states in 
strongly correlated magnets. 

RIXS cross section. In RIXS^ an x-ray excites an elec- 
tron from a core orbital to an empty state above the 
Fermi level. Rapidly the high-energy core-hole decays 
again. In this process the electron filling the core-hole 
emits an x-ray photon, which is then detected. In the 
short-lived intermediate state the motion of the valence 
electrons is affected by the presence of the core-hole po- 
tential. In case of indirect RIXS the perturbing core-hole 
interaction drives the creation of elementary excitations 
in the solid. Direct RIXS has a different scattering mech- 
anism: in that case the excited and de-excited electron 
simply involve different quantum states. The energy and 
momentum necessary to make excitations are lost by the 
scattered x-ray photons. By observing the energy loss u! 
and momentum loss q of the outgoing photons, one can 
thus determine the dispersion of the excited mode. 

The scattering amplitude Afi that takes an initial 
state |i) into a final state |/) is given by the Kramers- 

Heisenberg equation: Afi = J2n (/I ^ \^) ("I ^ N) / i^i ~ 
En + iT). The absorption of the incident photon and 
the emission of the outgoing photon are governed by the 
dipole operator l) cx A • p, with A the vector potential 
and p the momentum of the electron. The sum is over 
all intermediate states |n) with energies En and inverse 
lifetime F. In RIXS experiments, the energy of the pho- 
ton Ei is tuned to a resonant edge: Ei matches a set of 
En, which greatly enhances the scattering amplitude. 

In general the evaluation of the Kramers-Heisenberg 
equation is a complex problem because of the nature of 
the intermediate states. In many cases, however, the 
intermediate states are very short-lived, and one can em- 
ploy the lifetime h/T as a small parameter. This ap- 
proach is formalized i n the U ltra-short Core-hole Life- 
time (UCL) expansion.lSIinilll xhe Kramers-Heisenberg 
relation is expanded in H/T, where H is the Hamilto- 
nian of the system including the core-hole. The result 
is Af, = ^ {f\ DJ2Zo{iH - E,)/irf D\i) . Instead of 
evaluating the intermediate states, to first order one now 
only needs to compute the matrix elements of H. 

Single- and bi-magnon scattering. The experimental 
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effort to measure magnetic excitations with RIXS is fo- 
cussed on compounds in the copper-oxide family'^ To 
understand how three-magnon scattering comes about in 
these materials, we first briefly summarize the single- and 
bi-magnon scattering mechanism. Bi-magnon scattering 
at the Cu K-edge RIXS starts with an electronic Is — > 4p 
transition, with an energy of around 8990 eV. In the in- 
termediate state bi-magnon excitations are created dy- 
namically because the Is core-hole locally modifies the 
superexchange bond^^^. Single-magnon excitations are 
forbidden because the total spin is conserved in the scat- 
tering processes. For undoped cuprates the magnetic 
Heisenberg Hamiltonian in the presence of a core-hole 

where we restrict ourselves for simplicity to nearest 
neighbor exchange between spins Si at sites i, 5 denotes 
the neighbors and h\ is the core-hole creation operator. 
The first term is the spin 1/2 Heisenberg model with ex- 
change interaction J and the second term corresponds 
to the perturbation of the magnetic bond strength due 
to the presence of a core-hole. In first order of the UCL 
expansion one immediately obtains the two-magnon scat- 
tering amplitude A^^^''(q) cx {fl'^ie.^'^'^^h-Hhl |i),where 
Hi is the position of the copper ion absorbing the incom- 
ing photon and q the transferred momentum. It is im- 
portant to note that in the limit q — > selection ruleJ^i^ 
are such that the bi-magnon intensity rapidly vanishes 
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a q . 

In contrast, single-magnon scattering is allowed at 
the copper L- and M-edgefl^, where the initial photo- 
excitation is 2p —>■ 3d (L-edge, around 930 eV) and 
3p 3d (Af-edge, around 75 eV). The reason is that 
the spin of the p core-hole is not conserved in the inter- 
mediate states: it can flip via the spin-orbit coupling. 
The probability for spin-flip and non-flip are compara- 
ble since the core-hole spin-orbit coupling is very large, 
^ 20 eV at the L-edge. Therefore, when the core-hole is 
annihilated, a single-magnon can be created in the final 
statei^-i^. The spin-flip amplitude for a single copper ion 
is G(e, e') ~ (/I I'^i) i^jl ^ K) where rrij labels 
either the j = 1/2 or j = 3/2 spin-orbit split intermedi- 
ate states, depending on whether the L2- or La-edge is 
used, respectively. G{e, e') depends on the polarization 
vectors e and e' of the incoming and outgoing photons 
respectively. The resulting single-magnon scattering am- 
plitude is = _ (/I e^^-^^{St + S') \^). In 

the q — >■ limit, not only the single-magnon energy van- 
ishes, but also its spectral weight. Note that in addition 
to single-magnon scattering, at the L- and M-edges also 
bi-magnon scattering is allowed^ because the intermedi- 
ate state configuration is 3^^*^ and this singlet state com- 
pletely blocks any superexchange: J' = J in Eq. ([T]). But 
just as at the K-edge for q — >■ the bi-magnon intensity 
vanishes in leading order of the UCL also at the M- and 
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FIG. 1. RIXS three-magnon response function xai'i,^) for 
a two-dimensional Heisenberg model with nearest neighbor 
exchange. 



i-edges. Consequently at the Cu M-edge single- and bi- 
magnon scattering have negligible intensities: due to the 
low photon energy at the M-edge the transferred momen- 
tum |q| is restricted to the center 10% of the Brillouin 
zone of a typical cuprate. In second order of the UClJ^ 
a finite bi-magnon scattering amplitude appears at the 
A/-edge, but as at this edge typicalljisami j/p ^ g.l, this 
channel is very weak too. 

Three-magnon scattering mechanism. Single- and bi- 
magnon excitations can be created independently of each 
other, but in principle they can also be combined. Such a 
flipped spin in the vicinity of a bi-magnon excitation can 
form a tri-magnon bound state. We will show that the 
RIXS selection rule that forbids single- and bi-magnon 
scattering at q = does not apply for three-magnon 
scattering. Therefore three-magnon scattering becomes, 
in terms of the UCL expansion, a flrst order scattering 
channel. Consequently three-spin modes give a large con- 
tribution to the magnetic RIXS spectra of the high- Tc 
cuprates at the AI- and L-edges for small transferred mo- 
menta. 

The three-magnon excitation mechanism is very lo- 
cal and a direct combination of the single-magnon and 
a bi-magnon creation mechanism outlined above. It is 
caused by the 3(i^° intermediate state being screened by 
a bi-magnon excitation and then decaying via the spin- 
flip channel, which produces the third magnon. Via the 
spin-flip decay channel this process contributes to the 
flrst order of the UCL expansion, with an amplitude 

Af^ = (/I e'''-^'(^+ + S-)h,Hhl \^). Using 

the Hamiltonian in Eq. ([T]) and integrating out the core- 
hole degree of freedom yields the three-magnon scattering 
amplitude 

4? = ^^W^^^'''""" (5++5r)SrS.+. \^) . (2) 
This scattering operator can be written in terms of 
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FIG. 2. Calculated three-magnon L-edge RIXS spectra for 
the case that the [0, 7r]-direction is perpendicular to the scat- 
tering plane with scattering angles 29 — 90° and 130° respec- 
tively. In 90° geometry, not all momenta are experimentally 
accessible. 
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FIG. 3. Energy-integrated response functions Wi(q) for 
single-, two- and three-magnon excitations {i = 1, 2, 3, respec- 
tively). The ratio of the single- to three-magnon RIXS spec- 
tral weight is equal to (T / J)^Wi{q}/W:i{q}, see main text. 



magnon operators by expressing the spin operators in 
Holstein-PrimakofF bosons and then Fourier and Bogli- 
ubov transforming it. With conventions as in Ref. [Tl] 
the three-magnon scattering amplitude then becomes 

(3) JzG(e,e') 

k,p 



-7k+p+q (UpUk+p+q + VpV]^+p+q)] 
(/l4"i.«lk-p-qN)- 



(3) 



The cross section at zero temperature is obtained by 
summing over all final states with three magnons: / oc 
{Afil^ S{uj — ujfi) where w/i is the energy of the 
three-magnon final state /. We rewrite this as / oc 

— p4jY XsIQ;^); where X3(q, w) is the three-magnon 
response function. 

Three-magnon spectral shape. Fig. [T] shows the three- 
magnon response function X3(q, w) evaluated without 
magnon-magnon interactions. The three-magnon spec- 
tral weight is strongly peaked around uj = 5.2J and at 
q = 0. The reason for this is two-fold. First, the three- 
magnon DOS, which is a convolution of many combina- 
tions of three magnon states whose total momentum adds 
up to q, peaks close to its maximum energy of 3 x 2 J. It 
is in addition remarkable that the three-magnon response 
shows a strong q dependence, whereas the three-magnon 
DOS is quite featureless as the convolution tends to aver- 
age out dispersions. The strong q dependence is caused 
by the RIXS scattering matrix elements. Whereas these 
same matrix elements forbid at q = single- and bi- 
magnon scattering in first order, they act in the opposite 
way for three-magnon scattering, strongly enhancing its 
spectral weight at q = 0. 

A proper computation of for instance the tri-magnon 
boundstate in a two-dimensional antiferromagnet would 
require on top of this calculation the inclusion of three- 
spin magnon-magnon scatterin^i^Hi^l^ which is beyond 
our present scope. Considering magnon-magnon inter- 
actions in general terms, however, one expects these to 



reduce the three-magnon excit ation en ergy, in analogy to 
two-magnon Raman scatterin^iSIIIIlSI, xhe range of this 
energy reduction can be estimated by considering a Neel 
state of Ising spins: if three flipped spins form a string, 
the energy cost for such an excitation is 4J, whereas one 
isolated spin-flip and two neighboring flipped spins cost 
5J, three isolated ones 6J. Kinetics of course opposes 
three-magnon binding but in any case one expects a red- 
shift of tri-magnon spectral weight up to 2 J. 

Three-magnon polarization dependence. To unambigu- 
ously identify three-magnon spectral features in cuprates, 
knowledge of their theoretical polarization and scattering 
angle dependence is paramount. In evaluating the polar- 
ization factor G(e, e'), one is aided by the fact that in 
cuprates the magnetic moments are oriented in the Cu02 
planefPlSl]. Writing G(e,e') = e'^G^jC^ in Euclidian co- 
ordinates, the polarization factor for the Cu ^2,3- and 
Af2,3-edges simplifies to 



(4) 



If in an experiment the particular polarization of the 
outgoing photon is not detected, the cross section is 
summed over the two outgoing polarizations e'. When 
the [0, tt] -direction (in the Cu02 plane) is perpendicular 
to the scattering plane, and the incoming photons are tt- 
polarized, one obtainai^ G{e, e')^ cx cos^ (f>, where (j) 
is the angle of the incident beam with the normal to the 
Cu02 planes. This holds for all scattering angles 26. 

For the M-edge, changing the polarization vectors 
hardly affects q and therefore the spectral line shape 
is independent of the experimental geometry. It only 
affects the overall intensity, through G{e,e')'^. The Cu 
L-edge has a much higher energy: there is enough mo- 
mentum available to probe approximately the entire BZ. 
Fig. [2] shows the three-magnon RIXS spectra for the case 
that the [0, 7r]-direction is perpendicular to the scattering 
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plane with scattering angles 26 = 90° and 130° respec- 
tively. The asymmetry between +q and —q comes from 
the cos^ (j) polarization factor. 

Comparing one-, two- and three-magnon spectral 
weights. Because both excitations have the same po- 
larization factor, the ratio of three-magnon to single- 
magnon weight is independent of the polarization fac- 
tor G(e, e') within the UCL expansion. The RIXS in- 
tensity ratio is thus proportional to the ratio of the 
energy-integrated response functions of the two excita- 
tions, where the single-magnon response is given bjff^ 
Xi(q, u}) fa 0.718(Mq — Vq)'^S{uj — ujq). The numerical fac- 
tor ensures the sum rule (l/N) X]q/ duj xi(q, w) = 1 is 
obeyed, a physical property of the Heisenberg model that 
needs to be enforced within linear spin wave theory. The 
ratio of the RIXS single-magnon to the three-magnon 
spectral weight is equal to (T / J)^Wi{q) /W3{q) , where 
Wi{q) = J duj Xi(q, w). For typical cuprates at the Cu 
M-edge J/r « 0.1 and at the L-edgeP^ 0.2-0.4. Fig.[| 
shows the Wi{q) of single-, two- and three-magnon exci- 
tations. W^3(q) decreases by 8% from (0,0) to (tt, 0), and 
by 30% from (0,0) to (7r,7r). The bimagnon RIXS spec- 
tral weight scales with the polarization factor for non- 
flip scattering,'i^ so a detailed comparison of it to the 



single- and three-magnon intensity involves the experi- 
mental geometry of choice. In any case we can conclude 
from Fig. [3] that at q = and in a significant portion 
of the BZ around it the three-magnon spectral weight 
dominates over the other ones. 

Conclusions. We have theoretically demonstrated the 
importance of the three-magnon scattering channel in 
RIXS at the Cu L- and M-edge of the high-Tj, cuprates. 
For small transferred momenta, q — >■ 0, three-magnon 
scattering dominates the magnetic RIXS spectrum. This 
scattering channel gives direct access to the dispersion 
of for instance tri-magnon boundstates in gapped low- 
dimensional quantum magnets. Experimentally there is 
evidence for tri-magnon spectral weight in Cu M-edge 
RIXS data^on CaCuOs showing a single peak feature at 
an energy loss of 450 meV, consistent with our theory. 
This consistency can further be tested by measuring the 
polarization dependence and compare it to the theoreti- 
cal one presented here. Also in very recenlPS Cu L-edge 
data on La2Cu04 the three-magnon feature appears to 
be present around 400 meV for q = 0. These results 
illustrate that RIXS has an unprecedented capacity to 
probe multi-spin responses in strongly correlated mag- 
nets - both in principle and in practice. 
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